Discontinuous Percolation Transitions in Epidemic Processes, Surface Depinning in 
Random Media and Hamiltonian Random Graphs 



5h t 



Golnoosh Bizhani, 1 Maya Paczuski, 1 and Peter Grassberger 1 

1 Complexity Science Group, University of Calgary, Calgary T2N 1N4, Canada 

(Dated: April 9, 2012) 

Discontinuous percolation transitions and the associated tricritical points are manifest in a wide 
range of both equilibrium and non-equilibrium cooperative phenomena. To demonstrate this, we 
present and relate the continuous and first order behaviors in two different classes of models: The first 
are generalized epidemic processes (GEP) that describe in their spatially embedded version - either 
on or off a regular lattice - compact or fractal cluster growth in random media at zero temperature. 
A random graph version of these processes is mapped onto a model previously proposed for complex 
social contagion. We compute detailed phase diagrams and compare our numerical results at the 
tricritical point in d = 3 with field theory predictions of Janssen et al. [Phys. Rev. E 70, 026114 
(2004)]. The second class consists of exponential ("Hamiltonian", i.e. formally equilibrium) random 
graph models and includes the Strauss and the 2-star model, where 'chemical potentials' control 
the densities of links, triangles or 2-stars. When the chemical potentials in either graph model are 
O(logiV), the percolation transition can coincide with a first order phase transition in the density 
of links, making the former also discontinuous. Hysteresis loops can then be of mixed order, with 
second order behavior for decreasing link fugacity, and a jump (first order) when it increases. 
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I. INTRODUCTION 

Percolation describes the sudden appearance of 
system-wide connectivity arising from microscopic pro- 
cesses. It is a classic example [l| of a continuous ("sec- 
ond order") phase transition. Interest in systems where 
the percolation transition is discontinuous ( "first order" ) 
was sparked recently by claims for this in Achlioptas pro- 
cesses 0. Although these transitions were later shown 
to be continuous J-Q - albeit with unusual finite size 
scaling behavior [5|, the fact that percolation transitions 
could be discontinuous was claimed to be novel or even 
revolutionary @, 0, 0413|- After this, discontinuous per- 
colation transitions were observed in interdependent net- 
works [2514271 , in models inspired by but not using the 
Achlioptas trick p~5l4l8| | , and in a hierarchical lattice [2|| . 

One purpose of this work is to point out that discon- 
tinuous percolation transitions are not surprising and, in- 
deed, are common to a variety of (e.g. social or physical) 
cooperative phenomena. The existence of such transi- 
tions, together with an associated tricritical point, was 
proposed 25 years ago [2!| in the context of directed per- 
colation. This was verified in the seminal field theoretic 
work of Janssen et al. [30l | who introduced the gener- 
alized epidemic process (GEP) 0. In this scenario, the 
continuous transition is just ordinary percolation (OP), 



1 Janssen et al. called this the 'generalized general epidemic pro- 
cess' (GGEP), in order to distinguish it from the 'general epi- 
demic process' [D. Mollison, J. Royal Statist. Soc. B 39, 283 
(1977)]. The latter is now usually called an SIR (susceptible- 
infected-removed) epidemic, while a 'simple epidemic' in the no- 
tation of Mollison is now called SIS. We feel thus free to use the 
simpler acronym GEP for the generalized process. 



whilst the discontinuous one is the depinning transition of 
driven surfaces in random media at zero temperature [Ur- 
sa- 

Although the latter is continuous from the point of 
view of surface properties, it is discontinuous so far as 
the percolation order parameter is concerned 0. Indeed, 
the density of 'wetted' sites in the presence of a driven 
interface jumps discontinuously from zero to a finite value 
at depinning. 

A closely related line of papers finding discontinuous 
percolation transitions started independently in a social 
science context f39j - l4]| and addresses complex contagion 
or epidemics in random networks. It turns out that the 
model of [4(| is basically the random graph version of the 
GEP, as we explain in detail below. Our unified formu- 
lation based on GEP, that includes both social contagion 
and interface depinning, simplifies the description of both 
and isolates relevant variables that can affect the actual 
outcome in terms of potentially measurable observables. 

We present numerical simulation results for the GEP, 
including the time dependence of the number of growth 
sites in three dimensions. At the first order (= depinning) 
transition line, activity decays stretched exponen- 
tial in time whilst it behaves as a power-law both at the 
OP transition and at the tricritical point which separates 
rough from fractal growth. We give (tri-)critical expo- 
nents and compare them to theoretical predictions [301 ] . 

The main property that leads to first order transi- 
tions in the models we consider is cooperativity (or 'syn- 
ergy') in establishing links. This cooperativity can be 



2 Note that the co-appearance of first order bulk transitions with 
continuous surface transitions was discussed already in ISij . 
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implemented technically in different ways. We do this 
via stochastic dynamics as in [3(| H^-fSI for the GEP, 
and also via Gibbs-Boltzmann equilibrium distributions 
in Hamiltonian (or "exponential") ensembles, which have 
been used extensively to model social networks. 

Indeed, we also find discontinuous percolation tran- 
sitions in two exponential random graph models : the 
Strauss [H51 and the 2 " star model Q. They are both 
formulated in terms of a partition function and are gen- 
eralizations of the standard Erdos-Renyi (ER) random 
graph [45J . The Hamiltonians are bilinear with a control 
parameter (0) conjugate to the number of links and an- 
other control parameter conjugate to either the number 
of triangles or the number of "2-stars." When all con- 
trol parameters are 0(logN) (where N is the number 
of nodes in the graph) , the percolation transition can be 
either continuous or discontinuous, with hysteresis loops 
typical of first-order transitions. But for certain parame- 
ter regimes unusual hysteresis loops occur, where the per- 
colation order parameter exhibits second-order (singular 
but continuous) behavior for decreasing 9, but jumps dis- 
continuously for increasing 9. Such "mixed-order" hys- 
teresis loops have not to our knowledge been seen before. 

It is well known that the observation of continuity of 
a phase transition depends not only on the choice of the 
order parameter, but also on the choice of the control 
parameter. Take e.g. the standard example of a liquid- 
gas transition. If the temperature of water is increased 
at constant pressure, then the density and the free en- 
ergy jump discontinuously at the boiling temperature. 
If, however, the volume is kept fixed, no such jump is ob- 
served. Instead, as temperature is increased, a larger and 
larger fraction of the sample turns into vapor, but this 
happens in a completely continuous way. The standard 
assumption in thermostatics is that the order parameter 
is a density or inverse density (e.g. specific volume), and 
the control parameter is its conjugate (e.g. pressure). 
But in percolation, the standard choice of order param- 
eter is the fraction S max /N of sites belonging to the gi- 
ant cluster, while the control parameter is usually also 
a density - the density of occupied sites (bonds) in site 
(bond) percolation. Although this choice is legitimate, it 
can obscure the notion of first vs. second order transi- 
tions, since other choices more in line with thermostatics 
can lead to different conclusions. This might explain why 
previous works on first order percolation transitions were 
not recognized as such in the recent literature. 



II. THE GENERALIZED EPIDEMIC MODEL: 
COMPLEX CONTAGION TREATED AS A 
STOCHASTIC PROCESS 

Although the epidemic model of Ref. [30j| is formulated 
as a continuum field theory, the situation becomes more 
clear on a lattice. Consider a process where the probabil- 
ity of a given site becoming infected (or invaded) by one 
of its neighbors depends on the number of previous at- 



tempts by other neighbors. Once a site is infected, it tries 
once to infect every one of its not yet infected neighbors. 
Denote by pk the probability that an infection succeeds, 
if the attacked site has already fended off k previous at- 
tacks. If every attack increases the strength of the de- 
fender, pk decreases with fc, otherwise (if it weakens it), 
Pk increases. Site percolation is described by po > and 
Pk = for k > 1: If the first attack does not succeed, all 
later attempts are futile. Bond percolation is described 
by Pk = P for all k > 0. Ordinary (second order) perco- 
lation is observed whenever pk decreases with k (see also 
[4^|). but the transition switches to first order when pk 
increases sufficiently fast. In that case, infected clusters 
fill in most holes and bays, while protrusions are avoided 
- thereby making the clusters compact with rough but 
non-fractal surfaces. (The same effect is caused by high 
surface tension compared to disorder at the cluster-void 
interface in random media). Detailed predictions for the 
tricritical behavior in terms of an e = 5 — d expansion 
were given in (30j . 

A. The GEP on Random Graphs 

Although Dodds et al. [i(J EH assume a somewhat 
more complex mechanism of infection, their basic model 
can be mapped onto a sparse random graph model where 
each node with n neighbors in the the giant cluster is 
itself in the giant cluster with probability q n . This is 
precisely the mean field (random graph) version of the 
above model, if q\ = p and q n+1 = q n + (1 - q n )p n El- 
Due to the absence of short loops in this case, the condi- 
tion for tricriticality (transition between classes I and II 
in HO, El|]) simplifies to 

Q2=2qi, (1) 

with no restriction on any q n with n > 3. Since the 
derivation of this in [4(| EH is somewhat involved and 
obscures the relationship to the GEP as defined in [30j . 
a simple proof of Eq.([T]) is given in the supplementary 
material. 



B. The GEP on Regular Lattices: Tricritical 
Behavior and Rough Pinned Surfaces 

In the present work we studied in detail the case where 
Pk = p is the same for all k > 0, while po is different. 
Phase diagrams for simple (hyper-) cubic lattices and for 



3 Strictly spoken this is only true for static (i.e., percolation in 
the narrow sense) aspects at threshold. Dynamics and the be- 
havior above the threshold are different in both models, since the 
memory about previous contacts with infected neighbors is short 
lived in [4(1 l4l| , while it is long lived (previous attacks are never 
forgotten) in the model of [30(1 and in the simulations reported 
in the next subsection. 
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random regular graphs arc shown in Fig. 1. To the left 
of the curves, no infinite clusters exist, while such clus- 
ters do exist to their right. Since percolation thresholds 
on lattices scale as p c ~ l/(2d — 1) for large d, we used 
(z — l)qi and [z — l)y2 as coordinates in Fig. 1, where 
z is the coordination number. All curves start at site 
percolation (p = 0, 51 = 52 = Po)j since we do not con- 
sider here antagonistics effects (i.e., two attacks together 
cannot have less success than a single one). Tricritical 
points are marked by circles. There is no tricricital tran- 
sition in d = 2 [37]], i.e. isotropic rough 1-d surfaces are 
always fractal. For large dimensions, the lattice results 
converge to those for regular random graphs with degree 
z = 2d, as short loops become less and less important 
with increasing d. Preliminary results indicate that sur- 
face properties in the first-order regimes (i.e. below the 
tricritical points, but for qi > 0) are not in the same uni- 
versality class of surfaces without overhangs, which is the 
accepted theory for pinned rough surfaces |47H49| . while 
the situation is less clear in the case qi = 0, i.e. if at 
least two infected neighbors are needed for a site to be- 
come infected (in the latter case, epidemics can neither 
spread from single sites nor from (1, 0, . . .) surfaces, but 
they can spread from (1, 1,1...) surfaces). 

For the simple cubic lattice, simulations show that the 
tricritical point for this model is at at po = 0.111(2) and 
Pk = 0.464(8) for k > 0. Results of such simulations for 
epidemics starting from a single infected site arc shown in 
Fig. [21 where N(t) is the number of sites newly infected 
at time t. For OP, N(t) increases as a power law f with 
?/ pa 0.35 [lj], it decreases at the tricritical point as N(t) ~ 
t Vs with rj s = —0.70(1). This is in stark contrast to the 
prediction rj s ss 0.05 of [3(J. The tricritical point and OP 
are the only cases where N(t) shows a power law. For 
critical percolation with po > 0.111 the behavior crosses 
over to the OP scaling, while for po < 0.111 the data are 
compatible with a stretched exponential at the transition 
line (see lowest three curves in Fig. Analogous plots 
for the probability P(t) that the epidemic survives at 
least t time steps and for its average squared radius R 2 (t) 
are given in the supplementary material. They give S s = 
1.49(2) and z s = 1.205(4), where 6 S and z s are defined 
via P(t) ~ t~ s ' and R 2 (t) ~ t z ° . The predictions of (U 
arc S s ss 0.87 and z s ss 1.06. Again the agreement is far 
from perfect, although the changes from the OP critical 
exponents are in the right directions. More details are 
given in the supplementary material and in (50j . For the 
completely analogous case of directed percolation (SIS 
epidemics), see [5ll453| . 



III. COOPERATIVE PERCOLATION IN 
HAMILTONIAN RANDOM GRAPH MODELS 



The above discussion suggests that cooperativity in fi- 
nite temperature equilibrium systems may also lead to 
discontinuous percolation transitions. Indeed, the mean 
field percolation transition corresponds to the emergence 




FIG. 1. (Color online) Phase diagrams for the generalized 
epidemic process where infection of a site succeeds with prob- 
ability po in the first encounter with an infected neighbor, 
while it succeeds with chance p in all later encounters. Fol- 
lowing [4(J, we use instead of pk the probabilities qi — po 
and q-2 — po + (1 — Pa)p that that an infection has occurred 
after two encounters. Since percolation thresholds on regular 
graphs are roughly oc l/(z — 1) for large degree z, we use ac- 
tually [z — V)q n for the two axes. The curves labeled by d — 2 
to d = 6 correspond to d— dimensional hypercubic lattices, 
where z = 2d, while the curve labeled "z=15" is for random 
(i.e. locally loopless) networks. All curves start at the site 
percolation point qi — qi (small red circles). The bond perco- 
lation points (green crosses) are at 52 = (2 — 51)51. For large 
d they approach the tricritical points (big black circles) that 
converge as ((fc — 1)51 , (k — 1)52) —¥ (1, 2) for large k. The per- 
colation transitions are first order below the tricritical points, 
and second order above. For d = 2 there is no tricritical point, 
i.e. clusters and their surfaces are always fractal. In the first 
order regime (dashed curves), surfaces develop strong over- 
hangs and seem not to be described by models where these 
overhangs are neglected, except possibly at qi = (magenta 
squares), which seem to represent a different fixed point for 
d = 3. 



of the giant component in ER random graphs, where links 
appear independently with probabilit y p 1451 . The latter 
is the simplest "exponential model" [54j-|56|. In this ap- 
proach one considers graphs G with N nodes, were the 
probability for a given graph is defined by the Boltzmann- 
Gibbs equilibrium formula 



P(G; 6i,6 2 , 
withZ = £ G e- ff ( G - e 



-H(G;t 



(2) 



Here H is the Hamiltonian, 
{61, 62, ■ ■ ■} represents a set of control parameters, and we 
have set f3 = 1/kT = 1. More precisely, we assume that 
H is a sum of bilinear terms, 



H(G;6 1 ,6 2 ,...) =J2 e MG) 



(3) 
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FIG. 2. (Color online) Time dependence of the activity N(t) 
at five different pairs (po,p) where pk — p for all k > 1, for 
the generalized epidemic process on the simple 3-d cubic lat- 
tice starting from a single infected site. The uppermost (red) 
curve is for po — p — 0.2488 . . . which is critical bond per- 
colation. The middle (green) curve is for the tricritical point 
(po,p) = (0.111,0.464) (i.e., (qi,q 2 ) = (0.111,0.523)). Within 
the resolution of the figure, both these lines are straight. The 
lowest three curves are near the first order (depinning) tran- 
sition (p ,p) = (0.07, 0.792318 + A), where A = (magenta), 
+0.00022 (dark blue) and -0.00022 (light blue). Statistical 
error bars are all smaller than the line widths. 



where each A a is an observable ( "statistic" ) of the graph, 
and 9 a is the associated chemical potential. Typically, 
each A a represents the total number of small subgraphs 
(links, triangles, p-stars, 4-cliques, ...) in the graph. 

Models of this typ e have been popular in mathemati- 
cal sociology [ill, [5a] , although they tend to be unrealis- 
tic. In many cases, such models reduce to equivalent ER 
graphs without clustering and are trivial, apart from the 
usual, non-trivial dependence of the observables A a on 
the control parameters 9 a [57l |. 

The Hamiltonian for the ER model is: 



H ER (G;0) = 9L(G) 



(4) 



where L(G) is the number of links in G and 9 = ln[(l — 
p)/p\- It exhibits a percolation transition at p = 1/N 
(when N — > 00) [45| . thus the critical value of 9 is 

c = lniV. (5) 

In the following, we study the 2-star model [H| with 

ff 2 - 8 tar(G;0, J) = 9L(G) - ^n 2 (G) . (6) 

Here ^(G) is the total number of "2-stars", i.e. of pairs 
of links attached to the same node. We also consider the 
Strauss model 143, W3f with 



H stiauss (G; 9, B) = 9L(G) - ^n A (G) , (7) 



where Ha(G) is the total number of distinct triangles, 
i.e. of loops of length 3. In terms of the degree sequence 
{k h i = l...N}, 



while ua(G) depends also on degree correlations. 

For a typical (non-sparse) graph L increases quadrat- 
ically with N, while both n 2 and n& ~ iV 3 . This is 
why J/N and B/N are used as control parameters in 
Eqs. (|6I7[) instead of J and B. The 2-star model has, for 
any 9 > 2, a first order transition in the density of links at 
J* (9) and strong hysteresis. The results of [44| , together 
with a standard Maxwell construction show that 



r(9) 



(9) 



The line of first order transitions terminates at the crit- 
ical point 9* = J* = 2. (We neglect here all terms that 
are 0(1/ N) relative to the leading ones). Similarly, for 
the Strauss model a first order transition in the link den- 
sity p = (L)/N occurs for any 9 > 0.81 0. This time 
it is more complicated to obtain the exact transition line 
B*(9), but one can show that (see Supplementary Mate- 
rial) 



B*(9) w W for 0>1, 



(10) 



with a critical point at p* = 2/3, B* = 27/8 = 3.375, and 
9* = 3/2 - ln(2) w 0.807. 

For both models, a giant component exists in both the 
high and low link density phases, whenever 9 = 0(1). 
Thus the density transition happens when they are al- 
ready percolating, as long as 9 is finite. In order to reach 
a percolation transition one has to take 9 ~ IniV to get 
a sparse graph. In this regime the above estimates for 
the density transitions are still valid. Moreover, when 
B < B* or J < J*, respectively, the second terms in 
the Hamiltonians (Eqs. (|6I7[1 ) have no influence on the 
percolation transition, for N — > 00. This is illustrated 
for the 2-star model in Fig. 2, where we show numerical 
results averaged over 50 hysteresis loops for a network 
with N — 2000 and J = 3.0. The density transition 
(monitored via the average degree) indeed appears to be 
first order with strong hysteresis, while the percolation 
transition (monitored via (^max) /N) shows no hysteresis 
and is exactly the same as for ordinary ER networks. 

When J > IniV (or B > hxN, respectively) this sce- 
nario breaks down because the true equilibrium state at 
the ER percolation threshold, 9 p , is a dense graph that 
consists of a single giant component. Although the equi- 
librium network percolates at 9 p , a hysteresis loop start- 
ing at 9 > J begins with a sparse non-percolating graph, 
and due to metastability the effect of J (or B) is not seen 
until one passes the ER percolation threshold - provided 
that it remains in the metastablc region. This scenario 
is illustrated in Fig. 2 for J = 9 p = In AT. Now the hys- 
teresis loop for (k) /N is quite wide. The hysteresis loop 
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FIG. 3. (Color online) Hysteresis loops for the 2-star model 
when holding the chemical potential J fixed at a value of 
0(1) and sweeping the control parameter 6 for links. The 
loop on the left (red) is for the normalized average degree, 
while the one on the right (green) is for the percolation or- 
der parameter (S ma x)/-/V. The vertical dashed (blue) line is 
the Maxwell prediction for the true density transition. The 
curve for (S max )/Af shows practically no hysteresis and agrees 
within error with the one for OP. The percolation threshold 
9 C = In N is indicated on the x-axis. The rounding of the 
green curve near 9 C is a finite size effect (iV = 2000). 
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FIG. 4. (Color online) Analogous to Fig. 3, but for J = C = 
InJV = 7.601. This time the percolation threshold coincides 
with the true (Maxwell) density transition. As a consequence, 
also the curve for {Smax)/N shows hysteresis, with a contin- 
uous lower branch and a discontinuous upper branch. 



for {S max )/N shows the ordinary ER percolation shape 
on its lower branch, while it follows the discontinuous 
behavior of (k) /N on the upper (return) branch. To our 
knowledge, such a mixed-order hysteresis loop has not 
been observed before. 



IV. CONCLUSIONS 

As we already mentioned and as was pointed out re- 
peatedly before [H,0,[53], the 2-star and Strauss models 



are not realistic for real world applications. Accordingly, 
our demonstration that they exhibit both second and first 
order percolation transitions should be considered only as 
a proof that this phenomenon exists in equilibrium. More 
interesting examples can be easily suggested. A class 
of models that come into mind are random graphs (i.e. 
mean field type) where not only the number of nodes but 
also the number of links is set by hard constraints ( "mi- 
crocanonical models" ) . A model with two control param- 
eters (B and J) studied in [f38| exhibits zoo of metastable 
states, but it seems that most of these states are still too 
extreme to be physical. A further step in this direction 
could be to fix not only the average degree, but to fix the 
entire degree distribution, either by soft [59| or by hard 
[6fj| constraints. 

A more realistic class of Hamiltonian models with first 
order phase transitions could be spatially embedded (e.g. 
finite dimensional lattice) systems. Such models have 
not been studied much in the social science literature, al- 
though it is well-known that, for instance, spatial struc- 
ture is essential to maintain diversity in ecosystems. We 
believe that such models might provide a suitable mix- 
ture of structure and randomness to reveal important 
features of real, complex networks, including presumably 
percolation transitions of both continuous and discontin- 
uous type. 

In summary, we have shown that various percola- 
tion models can be naturally generalized such that they 
switch from ordinary, continuous behavior at the tran- 
sition point to discontinuous ("first order") behavior, as 
some parameter is varied. This parameter usually is a 
measure of coopcrativity in linking or "infecting" sites 
(the probability for sites to get linked is increased by 
other links already present), such that the percolation 
transition is more abrupt when cooperativity is high. 
We present a unified treatment including examples that 
range from social dynamics to condensed matter physics, 
and we also show that analogous phenomena occur both 
in stochastic dynamics out of equilibrium as well as in a 
Gibbs-Boltzmann equilibrium framework. We also sim- 
plify the dynamical description of the social contagion 
process introduced by Dodds and Watts [loj], clarifying 
thereby its relation to percolation and to the generalized 
epidemic process defined in [3(|. 

In condensed matter physics, the first order percolation 
transition is just the depinning transition of driven inter- 
faces in disordered media at zero temperature. Treat- 
ing it also in our unified framework not only allows us 
to study in detail the tricritical point (where we found 
for d = 3 striking disagreement with theoretical pre- 
dictions), but also to numerically investigate more effi- 
ciently a model for pinned surfaces in which overhangs 
are fully included and hence the rotational symmetry of 
the growth process is not explicitly broken at scales much 
less than the system size. In this latter context, the most 
important (but so far only preliminary) result we find 
is that overhangs are indeed crucial for such surfaces, 
and that all existing theories for critically pinned rough 
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surfaces (which neglect overhangs and are based on a being relevant to the most interesting physical case of 
single- valued "height function") might be obsolete, not isotropic media. 
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DETAILS ON THE SIMULATION OF 
TRICRITICAL AND FIRST ORDER 
PERCOLATION ON 3-D LATTICES 



All simulations were done on simple cubic lattices, with 
synchronous (discrete time) updates. We followed the 
spreading of epidemics that started either with point 
seeds or with the seed consisting of an entire infected 
plane. In the former we used lattices of size up to 2048 3 
and checked that clusters never reached the boundary. 
For small values of po, when growth from a point seed 
has a very high chance to die out, we used PERM [l| 
to grow clusters even if their probability was as low as 
10 -300 . For simulations initiated from an entire infected 
plane (results of which are not shown here but are used, in 
addition, to better estimate (tri-)critical points) we used 
lattices of sizes up to 4096 2 x 2048 with helical lateral 
boundary conditions. In that case we also implemented 
multi-spin coding in order to use only 2 bits to store the 
status of any site. We also recycled memory in order to 
grow epidemics that spread far from the initial infected 
boundary, by overwriting older parts of the cluster that 
were no longer growing. This enlarges the effective lat- 
tice size to 4096 2 x L z with L z > 2048. The precise 
thickness that we must not overwrite depended of course 
on the actual roughness of the growing surface. It was 
always checked that the cluster could grow without im- 
properly interfering with some of its older parts. 

Further details will be given elsewhere They will 
concern statistics, the precise methods used to estimate 
exact tricritical properties, critical exponents for rough 
pinned surfaces, and the behavior in dimensions different 
from 3. Here we present just two more figures (Figs. [1] 
and [2]), similar to Fig. 1 of the main text, that show the 
survival probability P(t) and the average squared radius 
R 2 (t) of newly infected sites, measured from the starting 
point of the epidemic. They show pwer law behavior at 
the tricritical point (straight lines on log- log plots), with 
exponent values given in the main text. 

A detailed comparison of our results with the field the- 
oretic predictions of 0] will also be given in 0] . Here we 
just mention that our results for the critical exponents 
are qualitatively similar (the changes relative to the ex- 
ponents for OP go in the right directions), but the agree- 
ment if far from perfect. The biggest disagreement is for 
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FIG. 1. (Color online) Log-log plot analogous to Fig. 1 of 
the main text, but for P(t) instead of N(t). Here P(t) is 
the probability that an epidemic started which with a single 
infected site is still growing after t time steps. The decay fol- 
lows a power law both at OP and at the tricritical point (with 
exponent —1.49(2) for the latter), while it seems to follow a 
stretched exponential at the first order (i.e., rough surface de- 
pinning) transition point. Again the central one of the three 
lowest curves (for po = 0.07) corresponds to depinning. 



T] s , defined via N(t) ~ t Vs . While a positive value, 



6 TT 45 



(1) 



with e = 5 — d = 2, was predicted in 0], we found r) s = 
-0.702(10). 



THE EPIDEMIC THRESHOLD IN THE 
DODDS- WATTS SOCIAL CONTAGION MODEL 

Since the original derivation of the result qi = 2q\ by 
Dodds and Watts [H is somewhat cumbersome and in- 
volves more than a minimal set of assumptions, we give 
here a simpler derivation following the typical arguments 
for epidemic thresholds via consistency conditions 043- 
We start by recalling the condition for the threshold of 
the standard epidemic process that leads to ordinary per- 
colation, and we then modify a suitable reformulation so 
that it also applies to the more general case with different 
infection probabilities pk- We finally obtain the critical 
line and the tricritical point by straightforward algebra. 

Let us call S the probability that a node at one end of 
a randomly chosen link gets infected during an epidemic 
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FIG. 2. (Color online) Analogous to Figs. 1 and SI, but for 
R 2 (t) which is the average distance of active (i.e. newly in- 
fected) sites from the seed. In the first order regime the clus- 
ter grows very slowly, and R 2 (t) should be dominated by the 
motion of their center of mass. Since clusters can grow only 
into new (not previously infected) areas, this motion should 
be essentially a self-avoiding random walk. Our estimate for 
the critical exponent is consistent with this, although there 
are large finite time corrections that make R 2 (t) grow even 
less fast than t for very long times, in particular for p below 
the transition. At the tricritical point, R 2 (t) is a power law 
with exponent 1.205(4). 

process on a sparse random network with degree distri- 
bution Pf.. If the infection can pass through any link 
with probability p, then the locally treelike structure of 
the network results in the consistency condition Q 



l-S = z- 1 Y,kP k (l-pS) k - 1 , (2) 
fc=i 

where z = (k) = J^k h^k- We write this as 

oo 

F OP (S) = J2 kPk{(l -PS)"- 1 + (S - 1)} = (3) 



k=l 



(the subscripts stand for 'ordinary percolation' 
percolation threshold is then defined by 



F OJ >(0) = F OP (0)=0 



The 



(4) 



where 

SI 



dFop/dS. Straightforward calculations 



give 



(5) 
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We also notice that F£ P (0) > 0. 

In order to modify this to arbitrary infection proba- 
bilities p n for attacks following n previous attacks, we 
first rewrite Eq. (|3|) such that contributions from differ- 
ent numbers of infected neighbors are separated. In order 
to do this we write 1 — pS = (1 — S) + (1 —p)S, such that 
the first term is the probability that the considered node 



is not infected, while the second terms is the probability 
that it is infected, but it cannot infect its neighbor since 
the link cannot be passed. Similarly we write 



fe-i 



Ef*- 1 )^-*)' 



k — n — 1 



(6) 



such that each term in the sum corresponds to exactly n 
infected neighbors. The modification to the generalized 
process is now obvious: We just have to replace the power 
(1— p)™ by (1— Po)(l— Pi) • • ■ (1-Pn-i)- Alternatively, we 
can replace it by 1 — q n where q n is the probability that 
n attacks succeed in infecting the site. The formulations 
using p n and q n are fully equivalent. 

Making this modification in Eq. (j3J) results in 



fc-i 



k=l 



n=0 



(7) 



{[l-q n }S n (l-S) 



k — n—1 



+ (s-i)}, 



where GEP stands for 'generalized epidemic process'. As 
before, the condition for criticality is 

F GEP (0)=i^ EP (0)=0, i^ EP (0)>0, (8) 

while the tricritical point is given by 

F GEP (0) = F^ EP (0) = F£ EP (0) = , (9) 

Evaluating the derivatives is straightforward and gives 



92 = 2qi 



(10) 



and qi is given by Eq. ([5]). Notice that the location of the 
tricritical point does not depend on any q n with n > 2, 
and its existence does not put any constraints on them. 

In the first order regime, the threshold condition for 
an epidemic is 



F GEP (S)=0, 



(11) 



which depends non-trivially both on the degree distribu- 
tion and on all p n (or all q n ). For regular graphs with 
degree z (i.e., Pk = 5k, z) and p n = p for all n > 1 
it approaches in the limit z — > oo the linear relation 
(72 = 1/z+pi. For z = 15, numerical solution of Eq. (11) 
gives the line plotted in Fig. 1 of the main text. 



THE STRAUSS MODEL 

We start from Eqs. (5) and (6) of Ref. 0, which read 
in our notation 



1 



-Bq(l-2/N) + I 



(12) 
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FIG. 3. (Color online) Panels (a) - (d) illustrate the cases 
(a) - (d) for the possible solutions of F(p) — 0. Panel (e) 
illustrates the Maxwell construction. The function 6{p) is 
given by Eq. (|19p . and the dashed horizontal line is placed 
such that the two areas between it and the curve 9 = 8(p) are 
equal. 



and 



1 + {e B / N - l)p 



(e*-B«(l-3/*0 + 1)2 + ( e B/N _ !j p 



(13) 



Here p is the link density and q is defined in Q . In the 
limit N — > oo and B,9 -C N of interest to us, these 
simplify to 



1 



P = 



-Bq + l 



(14) 



• (b) three different solutions (2 stable + 1 unstable), 
corresponding to phase coexistence (Fig. [3})); 

• (c) one single plus one doubly degenerate solution, 
corresponding to the boundaries of the coexistence 
region (Figs.[3tl and |3]:2); and 

• (d) one triply degenerate solution, corresponding 
to the critical point (Fig.[3Jl). 

The critical point is thus obtained from F{p) = 
F'(p) = F"(p) = 0, leading to the values quoted in the 
main text. The boundaries of the bistable region are 
given paramctrically by 



B = 



1 



= ln(l/p - 1) + 



1 



(17) 



2(1 -j>)p2 ' " *' ' 2(1 -p) 

which give asymptotically for large 9 

B-(6)&6, B+«exp(2(9) (18) 

for the lower (upper) boundaries in a plot of B versus 9 
(see Fig. 2 of If). 

The actual transition curve B*(8), or rather 9*(B), is 
obtained from a Maxwell construction: For some given 
value of B > B* , we first obtain 9 as a function of p from 
Eq. (HD, 



9{p) = Bp 2 +ln(l/p- 1) 



(19) 



For any 9q in the coexistence region, the equation 9(p) = 
9q has three roots p\ < pi < p%. The transition point 
9*(B) is then given by 



PA 



dp \9{p)-9*(B)] =0. 



(20) 



in 



For large B, the curve of 9{p) versus p tends (except near 
the points p = and p = 1) to a parabola, which gives 
in this limit 9*(B) sa B/3 as stated in the main text. 



and 



g= ( e ^ + l)2 = P 2 - ( 15 ) 

Combining these gives an equation for p in terms of B 
and 9 which we can write as 

F(p) = e e - Bp2 +1-1=0. (16) 
P 

This equation can have four outcomes (see Fig. [3]): 

• (a) one simple solution, corresponding to one single 
phase (i.e. no phase coexistence; Fig. [3£i); 
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